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Abstract 

The matrix sign function is shown to provide a simple and direct method to derive some 
fundamental results in the theory of surface waves in anisotropic materials. It is used to 
establish a shortcut to the basic formulas of the Barnett-Lothe integral formalism and 
to obtain an explicit solution of the algebraic matrix Riccati equation for the surface 
impedance. The matrix sign function allows the Barnett-Lothe formalism to be readily 
generalized for the problem of finding the surface wave speed in a periodically inhomo- 
geneous half-space with material properties that are independent of depth. No partial 
wave solutions need to be found; the surface wave dispersion equation is formulated 
instead in terms of blocks of the matrix sign function of % times the Stroh matrix. 



1 Introduction 



The theory of surface waves in homogeneous anisotropic elastic half-spaces has enjoyed 
remarkable progress in the 1960's and 1970's. The general framework developed by 
Stroh [l| for solving static and dynamic elasticity problems hasproved to be very fruit- 
ful for study of surface waves. In a series of classic papers, e.g. [2, H|, Barnett and Lothe 
employed the framework of Stroh to develop an elegant integral matrix formalism that 
underpins existence and uniqueness considerations for surface waves and that allows one 
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to determine the surface wave speed without having to solve for any partial wave solu- 
tions. The Barnett-Lothe integral formalism was quickly realized to serve as a corner 
stone for the surface wave theory On this basis, Chadwick and Smith [4| provided a 
thorough exposition of the complete theory for surface waves in anisotropic elastic solids, 
summarizing the major developments up to that time, 1977. Later on, a significant con- 
tribution came from Alshits. With his co-workers, he has done much work on extending 
the formalism of Stroh, Barnett and Lothe to various problems of crystal acoustics, see 
the bibliography in this Special Issue. A full historical record and a broad overview of 
the surface wave related phenomena may be found in [1, 0] . 

The purpose of this paper is both to present a fresh perspective on the integral 
formalism for surface waves and also to provide new results, including a generalization 
to laterally periodic half-spaces. The central theme is the use of the matrix sign func- 
tion which allows a quick derivation and a clear interpretation of the integral formalism 
of 0| . The application of the matrix sign function in the context of the Stroh formu- 
lation of elastodynamics was apparently first noted by |7J in the course of calculation 
of impedance functions for a solid cylinder. We reconsider the classical surface wave 
problem in terms of the matrix sign function, showing in the process that it provides a 
natural solution procedure. For instance, it is known that the surface impedance matrix 
Z satisfies an algebraic Riccati equation [8[, but it has not been used to directly solve for 
Z. Here we give the first explicit solution of this Riccati equation for the impedance. An- 
other important attribute of the matrix sign function is that it allows the Barnett-Lothe 
formalism to be readily generalized to finding the surface wave speed in a periodically 
inhomogeneous half-space whose material properties are independent of depth (i.e. a 
2D laterally periodic half-space). For this case, the present results provide a procedure 
that circumvents the need for partial wave solutions. Instead it establishes the disper- 
sion equation in terms of the matrix sign function which can be evaluated by one of 
the optional methods, in particular in the integral form similar to the Barnett-Lothe 
representation of the homogeneous case. 

The outline of the paper is as follows. The surface wave problem is defined in 
§|2] in terms of the Stroh matrix N. The matrix sign function signiN is introduced 
and discussed in $3] where it is shown to supply a novel and possibly advantageous 
route to the Barnett-Lothe integral formalism. An explicit solution of the algebraic 
Riccati equation for the surface impedance matrix Z is derived in §HJ Application of the 
matrix sign function to formulating and solving the surface wave dispersion equation 
in a laterally periodic half-space is considered in £j5j with the numerical examples given 
for a bimaterial configuration. The Appendix highlights explicit links between the sign 
function and some related matrix functions. 

2 Background and problem definition 

The equations of equilibrium for time harmonic motion (with the common factor e~ lujt 
everywhere omitted but understood) are 



where p is mass density, are the elements of the elastic tensor C referred to an 
orthonormal coordinate system, Ui and are elements of the displacement u(x, t) and 
stress cr(x, t). We first consider a uniform elastic half-space n • x > 0, (|n| = 1) with 




pu 2 Ui, cTij = Cij k idiu k (i,j,k,l 
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constant density and elastic moduli. Solutions are sought in the form of surface waves 
propagating in the direction m parallel to the free surface (m ■ n = 0, |m| = 1): 



u \ _ / a(n-x) 
n<r J \ — ik\{ri ■ x] 



e ifcmx (2) 



The equations of equilibrium fll]) take the form of a differential equation for the 6-vector 
10, 

|| = ikNr,, with n{y) = ffiff) , (3) 

N = _ ( (nn)-\nm) (nn)- 1 \ ,^ 

Umn) (nn)" 1 (nm) — (mm) + pv 2 I (mr^inn)' 1 J ' 

where v = u/k, the 3x3 matrix (pq) has components (pq)jk = PiCijkiqi for arbitrary 
vectors p and q, and I is the identity matrix. The real-valued Stroh matrix N satisfies 
N T = KNK where T indicates transpose and K in block matrix form comprises zero 
blocks on the diagonal and identity blocks off the diagonal. Denote the eigenvalues 
and eigenvectors of zN by \ a and £ Q = (a Q , l a ) T (a = 1, 6), and introduce the matrix 
r = H^, £ 6 || with columns £ a . Assume in the following the normal situation where all 
A a are distinct. Then the above symmetry of N yields the orthogonality/completeness 
relations in the form [1, eqs. 2.8-2.12] 

r T KT = 1 TKT T = I, (5) 

where the normalization £^K£ Q = 1 has been adopted. Hereafter we use the same 
notation I for the identity matrix regardless of its algebraic dimension. 

Throughout this paper we restrict our interest to the subsonic surface waves and 
thus assume that v is less than the so-called limiting wave speed v, see [1, 0, 0]. This 
implies that i\ a are all nonzero and in pairs of complex conjugates (denoted below by 
*), so the set of eigenvalues and eigenvectors of zN can be split into a pair of triplets as 

A q = -A; +3 (ReA Q < 0) C = ta+3, " = 1,2, 3. (6) 

These two triplets are commonly referred to as physical (a = 1,2,3) and nonphysical 
(a = 4, 5, 6) since they define partial modes that, respectively, decay or grow with 
increasing n • x. The eigenvector matrix T partitioned according to (j6j) is written in the 
block form as 

Ai A 2 \ Ai = ||ai,a 2 ,a 3 || , A 2 = ||a 4 , a 5 , a 6 || = A*, 

with (7) 
L i L 2/ Li = ||li,l 2 ,l 3 || , L 2 = ||1 4 ,1 5 ,1 6 || = L*, 

where the blocks A 1; Li and A 2 , L 2 describe the physical and nonphysical (decaying 
and growing) wave solutions, respectively. From ([6]) T* = TK and so the orthogonal- 
ity/completeness relations (jSJ) may be cast in the subsonic domain to the form 

r+KT = K TKT + = K, (8) 

where + denotes the Hermitian transpose. Note that the relations (jHJ) admit interpreta- 
tion in terms of the energy flux into the depth. 
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The surface wave solution comprises the decaying solutions only and therefore must 
have the form 

a(y) = Ax e yX v, l(y) = L, e yX v with A = diag(A 1; A 2 , A 3 ), (9) 

where v is some fixed vector. The surface wave problem for the homogeneous medium 
is posed as finding for a given wavenumber k G M + the surface wave speed v — v s £ R 
for which the surface traction 1(0) vanishes, and hence v is a null vector of Lj (although 
this is not a fruitful avenue to follow, which is the whole point of the Barnett-Lothe 
solution procedure based on the integral representation). 

A variety of related notations have been used for the surface wave problem. We 
generally follow Barnett and Lothe 0] where the notation is based upon that of [jjj]. 
Barnett and Lothe 0] also provide comparisons of their notation with that used by 
Chadwick and Smith [4], which is closer to that of Ting [5]. The slight notational 
differences are related to the choice of the 6— vector in ([3]), and amounts to different 
signs for the diagonal or off-diagonal elements of the 6x6 matrix analogous to N. 



3 The matrix sign function signzN 



3.1 Definition 

The sign function of a matrix M is conveniently defined by analogy with the scalar 
definition as 10 



signM = M (M 2 ) 



-1/2 



(10) 



where the principle branch of the square root function with branch cut on the negative 



real axis is understood; z 



-2\l/2 



signz with signz = +1(— 1) if Rez > 0(< 0). As a 



result, the sign function of a matrix M with eigenvalues and eigenvectors denoted by X a 
and £ Q satisfies 



(signM)£ Q = ±£ Q for X a G C 



± 



'ir 



Note that signM is unchanged under M — > aWL+ibl, a £ IR + , 6g1, and it is undefined 
for eigenvalues lying on the imaginary axis (ReA a = 0). We also note for later use the 
property sign (UMU -1 ) = U( sign M)U _1 . Evaluation of the matrix sign function is 
possible with a variety of numerical methods 11] . the simplest being Newton iteration of 
M fc+1 = |(M fc + M fc 1 ), M = M, with M k -»■ signM in the limit as k ->• oo, although 
this can display convergence problems. Schur decomposition, which does not require 
matrix inversion, is very stable, and is readily available, e.g. [12]. The function signM 
also has integral representations, which we will use in order to shed fresh light on the 
integral formalism in the surface wave theory. See HI, [ll] for reviews of the matrix sign 
function. 

The following expression for the matrix sign function is based o n (flOl) combined 
with an integral representation for the matrix square root function 14], 111 ] 



signM = -M dt(t 2 l + M 2 y\ 
71 Jo 



'12) 
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Equation ([T21 may be converted into the following form using a change of integration 
variable 

sign M = (M e ) with M e = (cos 01 - isin0M) -1 (cos0M - z sin 01), (13) 

where 



7T _ _ 

denotes the average. Differentiation of M# yields 

dM 



d0 



i(M 2 e -l), (14) 



which implies (e J<?iMs ) = (e *^ M9 ) 1 = cos <pl + z sin <p(M.$) for any 0. This provides 
alternative identities for the matrix sign function, such as 

. / S in(0M e ) x Mg cos(0M e ) x 

signM = ( — — ) = ( ), 15 

v sm0 cos0 

where the value of <fi is arbitrary as long as the denominator, sin0 or cos</>, does not 
vanish. The special case of eq. (IT3l i corresponds to the limit as <f) — > in (|T5l) . 



3.2 signzN and the Barnett-Lothe matrices 

Let us apply the above definition and properties of the matrix sign function to the case 
sign z'N where N is the Stroh matrix given in (j3J) and considered in the subsonic domain 
v < v. The eigenspectrum of N is assumed partitioned according to (J6]). Hence equation 
(fTTjl o taken with M = z'N reads 

(signz'N)£ a = -£ a , (sign z'N) £ Q+3 = £ a+3 , a = 1,2, 3, (16) 
Denote the blocks of sign z'N (= K (sign z'N) J K) as 

sign z'N = i (t qt^J with Q = Q T , B = B T , trS = 0, (17) 

where S, Q, B are all real since so is N. Note that sign z'N and hence S are traceless. 
Appearance of the Barnett-Lothe notations on the right-hand side will become clear 
in the course of the upcoming derivation. The involutory property of the matrix sign 

function, ( sign z'N) = I, implies the identities 

BS + S T B = 0, SQ + QS T = 0, S 2 + QB = -I. (18) 

Taking into account the spectral representation z'N = Tdiag(A, — A*)T _1 and the relation 
©2 yields the spectral decomposition of the matrix sign z'N in the form 



signzN = rdiag(-I,I)r T K= 2 2 _ 1 J T . T _. T 



A 2 L«2 — AjLf A 2 A 2 r - AiAf 



(!9) 

AiLf AiAf \ _ JAA ( A l 
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where we have noted the projector P_ on physical modes: P T = (see more in 

Appendix). From (JTZj) and (CHI . 

S = z(2AiLf - I), Q = 2iAiAf , B = 2iLiL^. (20) 

Assume that Ai and hence Q are invertible (we will return to this point later). 
Introduce the surface impedance Z = iLiA^ 1 . It is Hermitian due to (JHJ) 1 - Plugging 
([IT]) into (|T6l) i gives SAi + QL X = iAi which, with regard for the invertibility of Q, 
enables expressing Z through the blocks of signzN. Thus 

Z = ^LiA^ 1 (= Z+) ^ Z = -Q- x (I + zS) . (21) 

Note from (J2"T]) that, for any v, det S = (since Q _1 S is antisymmetric, see also (ITSl) -?) 
and that det (I + iS) is real (since Z = Z + ). 

General significance of signzN for surface waves is immediate when one considers 
the boundary condition which demands that a non-zero surface displacement exerts zero 
surface traction. Since a surface wave is composed of decaying modes only, this means 
there must exist a vector v^O such that 

' (I + iS)v = 
J V J " I Bv = 



8ign»N)(^)= n V =► r""""'' (22) 



Hence both detB and det (I + iS) must vanish at the surface wave speed v — v s . For 
the case in hand of a uniform material, these are not two independent conditions but 
one condition. Indeed, the traction-free boundary condition can be posed in the form 
det Li = 0, hence from (1201) and (j2T]) i determinants of B, of I + iS (= 2A 2 L^) and of Z 
vanish simultaneously. Note that the equation det (I + iS) = reduces to tr (S 2 ) = —2 
(since trS = and det S = 0). Also at v = v s , the common null vector v of I + iS and 
B can be cast as v = a + ib with linear independent real vectors a and b, hence B is 
rank one. To summarize, the subsonic surface wave speed v = v s can be determined 
from det Li = or else from any of the equivalent real dispersion equations (see (a, eqs. 
(12.10-l)-(12.10-4)]) 

det Z = tr (S 2 ) = -2 det B = (trB) 2 = tr(B 2 )) , (23) 

which all are expressed through the blocks of the matrix signiN. 

Equations (flB]) . ( 1201) . ( 12T1) . ( 1231) are basic relations of the Barnett-Lothe formalism 
with regard to surface wave theory. Interestingly, we have arrived at these relations 
by single means of the definition of sign function of i times the Stroh matrix N = 
KN T K, without specifying the method of this function evaluation and without explicitly 
attending to the fundamental elasticity tensor which underlies development of Barnett- 
Lothe formalism in @, d, 0]. It is instructive to highlight an explicit link between the 
two lines of derivation, see next. 

3.3 Linking sigruN to the fundamental elasticity tensor N# 

Following 0,0,0], introduce the so-called [3] fundamental elasticity tensor 



N, 



\rs]\s8\-Hirr] - kri Irsllss}- 1 



where \pq\ = (pq) — [m ■ p] [m ■ q] pv 2 l with p, q = r, s : 
r = cos #m + sin n, s = — sin 9m+ cos 6n. 
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The matrix N# is defined in the subsonic domain which means that v is small enough 
to guarantee existence of [ss] -1 for any 9. Note that at 9 = is equal to the Stroh 
matrix N defined in 

Lemma 1 The matrix sign function ofi times the Stroh matrix N is i times the angular 
average of the fundamental elasticity tensor N# 

S1 gn,N = ,(N e ) B = <M-H[Vh>. (25) 

Proof: Equation ( TT3]) with M = zN reads 

signiN = i(G) with G = (cos 91 + sin 0N) _1 (cos 9N — sin 91) . (26) 



Let us show that N# defined in (EMI) and G given by ff26]) 9 are one and the same. 
Substitution into G the representation [3, eq. (3.15)] 



N = (X- Yr 1 J(X + Y) where J 



I 
-I 



x / ( 27 ) 

nn) \ _ / [mm] 



X-Y= > ' " , X + Y- i _ 
\[mn) —I J y—{nm) — 1 

yields, using e je = cos 01 + sin^J, 

G = ( cos 0(X - Y) + sin 0J(X + Y)) -1 ( cos 0J(X + Y) - sin 0(X - Y)) 
= ( e ™X _ e - je Y) J (e je X + e- je Y) 

= (X - e- 2je Y) 1 J (X + e- 2je Y) . (28) 
Noting the identities 

X_e Y -{lrsj -ij' X + e Y "^-[ S r] -ij' (29) 
it follows that 

Combining f[3"Uj) with (J2"5]b completes the proof of (125]) . ■ 

Note that by the definition of G in ( |26l) the eigenvectors of N#(= G) are indepen- 
dent of 9 and that eq. (|14p with M# = zNg reads 

^-I-NJ. (3D 

Thus we have used the defining properties of sign function of the Stroh matrix N to 
obtain the basic relations f[T8"j) , fl2"0"j) . fl2"Tj) . fl2"3"|) of the Barnett-Lothe integral formalism, 
and then we made appear the fundamental elasticity matrix Ng but only in the context 
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of one of possible (integral) representations of the sign function. This is different from 
the original methodology of [HTH, 0| which follows a somewhat inverse line of derivation. 
It proceeds from Ng and establishes its properties such as fT5T]) which are then used to 
arrive at the relations ( 1181) . ( 1201) . ( l21i) . ( 123]) . One way or another, the core is that the 
surface wave speed v = v s can be defined through the matrix sign function which admits 
evaluation by different integral formulas, e.g. eqs. ( fl2|) . ( TTBT) . ( IBTj) and ( |69|) . or iteration 
schemes, or otherwise. Interestingly, some of these 'indirect' methods for calculating 
(N#) have been discussed in the literature, thou gh with no mention of the matrix sign 
function. For instance, Gundersen and Lothe 15| propose a scheme tantamount to 



Newton iteration, while Condat and Kirchner 161 describe a method based on continued 



fractions that appears to be similar to a stable algorithm for computing the matrix sign 



function 17 



In conclusion, a remark is in order concerning invertibility of Q in the subsonic 
domain v < v. By definition of the latter, the eigenvalues N for v < v cannot lie on the 
real axis, so G in (126]) -? is well-defined and hence, due to G = N#, the matrix [[ss]] in 
(|24|) is non-singular. Since [[ss]] is obviously positive definite at v — 0, it follows that 
[[ss]] is positive within the whole subsonic domain v < v (this of itself is an alternative 
definition of the subsonic domain). Using Q = — ([[ss]] -1 ) as given in (125]) shows that 
Q is indeed invertible in the subsonic domain which implies that A 12 are invertible and 
Z exists at any v < v. Note that we did not employ the fundamental elasticity tensor 
Ng in this proof. 



4 The surface impedance as defined by Riccati equa- 
tion 

The surface impedance matrix Z, which is introduced in the subsonic domain by any of 
the following relations 

\ a = -iZa a (a = 1,2, 3) <=> l(y) = -iZ&(y) <=> Z = iL 1 A^\ (32) 

plays a central role in the theory of surface waves. It is a crucial ingredient in surface- 
wave uniqueness and existence considerations [!, HI] and is of fundamental significance for 
energy and flux l|| • Direct use of the impedance has found wide application, including 



guided waves and scattering in multilayered structures, e.g. [19|, |20|, 12JJ, |22] , edge waves in 
anisotropic plates and cylindrical shells [HI, [HJ , and waves in cylindrically anisotropic 
solids [Tj. The surface impedance Z has been shown in §3.21 to satisfy the Barnett- 
Lothe formula ( 121]) which can be evaluated by integral representation or other optional 
expressions for signiN. Alternatively, the impedance can be defined via differential 
and algebraic Riccati equations. Thus, Biryukov [H, H developed a general impedance 
approach for surface waves in inhomogeneous half-spaces based on the differential Riccati 



equation, see also [26]. For the homogeneous half-space considered here, substitution of 
1 = — iZ& in eq. ([3]) combined with the fact that the impedance Z associated with the 
surface wave solution (jUJ) is a constant matrix independent of y leads to the algebraic 
Riccati equation [1, eq. 3.3.7] 

(Z — i(mn)^ {nn)~ l (Z + i(nm)) — (mm) + pv 2 I = 0. (33) 

This equation has a unique Hermitian solution for Z which is positive definite for v < v s 
2l\ . It is clear that the matrix Z of the Riccati equation (j33J) must be identical to Z 
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defined within the Barnett-Lothe formalism, see [27|; however, no explicit solution of 
the algebraic Riccati equation has been provided. Here we derive such a solution to ( 1331) 
using a method based on the matrix sign function (in fact it is for this class of matrix 
problems that the sign function was first proposed [14j). 

The main idea of this method is to rewrite the Riccati equation in a form where 
the matrix sign function reduces it to a linear equation for the unknown matrix: 

Lemma 2 For some matrix R the algebraic Riccati equation (1331 is equivalent to the 
following identity, 



(35) 



JZ OJ \ I J \iZ / 

Substitution from ([T7|) i into (!34l) yields 

S Q\ /-I I\ /-I i\ /-I R 
B S T IjZ o Uz oM I 



hence, znter aZia R = I — iS, and more importantly 

Z = -Q -1 (I + iS), B = (l + iS T ) Z = Z(I-iS) (36) 

where ( 13"6l) i is the same as ( l2Ti) and ([361)2,3 follow from combining identities ( fl6|) i and 
(ITS]) (note that I + iS T and I - zS are not invertible at the surface wave speed v = v s ). 
We have therefore constructed an explicit solution of the algebraic Riccati equation ( 1331) 
for the surface impedance tensor Z and shown that this solution coincides with the 
Barnett-Lothe expression (12 lj) . It remains to justify the claimed identity of Lemma EJ 
Proof of Lemma [2j Consider the 6x6 matrix equality 

•-ivt Y N i N A f- 1 A ( M i M \f-i iv 1 . . 

where 

M = Z" 1 N 3 , M 1 = -Z _1 (N 3 + iN? , Z), M 2 = Z _1 (N 3 - iZNi). (38) 

Simple algebra shows that the matrix Riccati equation (133!) is identical to the upper 
right block in (1371) while the other blocks are trivial identities. Referring to the matrices 
Mi and M 2 , we note that the eigensolutions X a and £ a = (a a , \ a ) T (a = 1, 2, 3) satisfy, 
from eq. (j3J), 

- N 3 a a + N[1 Q = i\ a \ a =► -N 3 Ai + Nf Li = iUX. (39) 

Replacing L x = — zZAi in (139^) yields Mi while M 2 follows similarly from the substitu- 
tion Ai = zZ _1 L 1; 

Mi = AiAA^ 1 , M+ = -LiAL^ 1 . (40) 

It is clear from these representations that Mi and M 2 have eigenvalues in C~ and C + , 
respectively. Define R as the solution of the following Sylvester equation 

I ±R\~V M i M \ fl |R\ /Mi 



RM 2 -M lR = 2M ^ .j ) ^< -\) = { M J, (41) 

then ( 13"4"|) follows from eq. (1371) and signdiag(Mi, M 2 ) = diag(— I, I) 
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5 Surface waves in a laterally periodic half-space 

5.1 General setup 

The equations of equilibrium (jTJ) in a spatially inhomogeneous half-space with p = p(x), 
Cijki = Ctjfci(x) may be written in Stroh-like format reminiscent of eqs. ((3]) and (J4j) as 
follows. Distinguish the coordinate normal to the surface, y = n • x, and the vector 
z G M? representing the orthogonal two-dimensional space. The governing equations 
then take the form 

^ = iAf/i (42) 

with 

/ u \ / T' l TZ + T' 1 \ 

z) = [ixur ) ' M = ~ [tzt- 1 ^ -q + P lu 2 i nr^ 1 ) ' (43) 

T = (/m), K = HV 2 , n), Q = {-iV z , -iV z ), 

where the definition of the 3x3 matrix (pq) = (p,q) has been expanded to include the 
operator with components (pq)jk = PiCijkiQi for vector operators p and q. 
Let us consider the material with parameters that are T-periodic in z: 

h(y, z + n x ai +n 2 a 2 ) = h(y,z) for h = p,c ijkh (44) 



where z G M 2 , n 1 ,n 2 G Z, and the linear independent translation vectors a 1; a 2 G 



d2 



define the irreducible unit cell T = Y^>=\ e IP> 1]) °f the periodic lattice. Let 



{e^} define an orthonormal base in R and denote 

a,- = Ae j; b, = (A- X ) T e j; T = {g : g = £^27171^, (m, n 2 ) G Z 2 }. (45) 
The material parameters can therefore be expressed as 

%,«) = y)e <g -*^(y,g) ^ %,g) = -L [ dze-^ z h(y,z) ioih = C,p. (46) 
g6 r 1 1 1 

In practical terms the matrices are limited to finite size by restricting the set of reciprocal 
lattice vectors to g G To where To C T comprises a finite number of elements, say N. 
We look for solutions in the Floquet form 

/i(y, *) = n(y,z) e*"- (47) 

for T-periodic r)(y,z). The surface wavenumber vector k = /cm resides in the first 
Brillouin zone of the reciprocal lattice T and is otherwise arbitrary. Employing a plane 
wave expansion 

r7(y,z) = ^ e**f?(l/,g), (48) 

eq. (14*2]) becomes an ordinary differential equation for the 6 N— vector rj(y) comprised of 
all r)(y,g), g G T , 
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Here P = P + , T = T + and R are 3N x 3iV matrices with 3x3 blocks associated with 
PWE wavenumbers g, g', given by 

(P) g g' = (* + g,k + g') (g - g') - u 2 p(y , g - g')I, 
(T) gg , = (n,n)(g-g'), (RW = (fc + ^n)(g-gO, 

where the definition of (pq)jk(g) is the natural extension of (pq)jk to include the Fourier 
transform, 

(pq)jk(g) =PiCijki(y,g)qi- (51) 

Note that T is negative definite and hence invertible. 

Equation ( |49p is valid regardless of whether the material properties depend on y 
or not. We consider here the case of a laterally periodic material whose properties are 
periodic along the surface and uniform in the depth direction, so that the system matrix 



N is constant. A method for treating the case of periodic N(y) is discussed in [28 



Surface impedance and related matrices via the matrix sign 
function 



5.2 



For constant N, the governing equation (1491) is analogous to that in the uniform elastic 
half-space and its solution is in a way analogous to ([9]), with two major differences. The 
first is that we are now dealing with large, formally infinite, vectors and matrices (for 

convenience, we continue to refer to N of 6N x 6N size). The second is that, in contrast 

to the real-valued Stroh matrix N, the system matrix N of ( 1491) is generally complex. 
Its symmetry is 

N = KNK, (52) 

where K is a 6N x 6N equivalent of the 6x6 matrix K with two zero and two identity 
blocks on and off the diagonal. Denote the eigenvalues and eigenvectors of zN by X a 
and £ Q As everywhere above, we restrict our attention to the subsonic domain where, 
by virtue of ( 1521) . the set of 6N eigenvalues of zN can be split in two halves as 



A, 



-A 



a+3N 



(ReA Q < 0) for a 



1,...,W, 



(53) 



which correspond to physical (decaying) and nonphysical (growing) modes. Adopting 
the same partitioning for the eigenvectors, denote 



6JV 




(54) 



It also follows from f[5"2"j) that the orthogonality /completeness relations in the subsonic 
domain hold in the form similar to (IE]), namely, 



r+KT = K <£> TKT 



K. 



(55) 



Note that, in contrast to the uniform case ([6]), £ a and £ a+3N are not complex conjugated, 

i.e. r* TK and so f l5l)|) has no equivalent that would be a generalization of Q. This 
is reminiscent of the Stroh formalism in cylindrical coordinates, see [7| and [29[ §3d]. 
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Now we can follow the line of derivation proposed in §3.21 Introduce the sign 
function associated with zN, namely, 

(signiN)^ a = -l a , (signzN)£ Q+3 = l a+3 , a = 1, ...,3N (56) 
and denote its blocks as 

signiN = z^| S\ with Q = Q+, B = B+ tr S = 0. (57) 

By definition ( 1561) . (signiN) 2 = I and so 

BS + S+B = 0, SQ + QS + = 0, S 2 + QB = -I. (58) 
Using (|55|) . the spectral decomposition is (cf. eq. (119]) ) 



sign »N = f diag(-I, I)Kf +K = °) + 2 (f J ? gg) . (59) 

Hence from (1571) and (159)) 

S = i(I - 2A 2 L+), Q = -2iA 2 A+ B = -2iL 2 L+. (60) 

Assuming invertibility of Q and hence of Ai and A 2 , introduce the surface impedance 
Z = iLqAj -1 . From ( )55|) Z is Hermitian and combining ( )56|) i with ( 157|) relates Z to the 
blocks of signiN. Thus 

Z = zL 1 A^ 1 (= Z+) ^ Z = -Q-^I + iS). (61) 

Note from (IBT]) that det (I + iS) is real (since Z and Q are Hermitian) and, unlike the 

case of a uniform medium (hence of real Q and S), the determinant of S is generally 
nonzero. . _ 

Finally, the evaluation of signzN is possible by different methods, including the 
integral representation analogous to ( 125|) . yielding 



S = -<T e - 1 R+>, Q = -<T- 1 ),^B = <T e+f -R e T ! 1 R+>, (62) 
T e = cos 2 6T + sin 2 6P - sin 9 cos 0(R + R + ) , 
R = cos 2 0R + sin 2 #R + + sin cos 0(T + P). 



where 



5.2.1 Dispersion equation 

Existence of a surface wave with a speed v — v s under the traction-free boundary implies 
that there must exist some non-zero vector v at v s such that 

• >i*(v\ (-A /(? + <S)v = 0, /A 2 L+v = 0, 

srgnzN ~ = ~ < ~ ^> < (63) 

\ Bv = 0, 1l 2 L+v = 0, 
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where (|57|) and f[6"Uj) have been used. The situation is reminiscent of eq. (I2"2"j) for a 
uniform half-space, certainly apart from the fact that eq. ( 1631) involves matrices of a 
large, formally infinite, size. Another particularity is that the equality det B = may 
signify occurrence of either a physical (if det Li = 0) or a nonphysical (if det L 2 = 0) 
surface wave solution. At the same time, the condition det(I-MS) (= det2A 2 L+) = is 
both necessary and sufficient for the physical wave in, specifically, the subsonic domain, 
where Q and Ai, A 2 are invertible (and so Z exists) as can be argued similarly as in 
§3.3. Interestingly, this is no longer so in the upper overlapping band gaps of the Floquet 
spectrum, where the condition fl63|) is recast in the form with a positive definite matrix 

((I + 2S) + (I + *S) + B + B)v = (<£> L 1 (A+A 2 + L+L 2 )L+v = 0), (64) 

which provides a single dispersion equation 

det ((I + iS) + (I + iS) +B+B) = (^ | detL^ 2 = 0). (65) 

Note that those considerations above which used spectral decomposition under the as- 
sumption of distinct eigenvalues can be reproduced in the invariant terms of the projector 
matrix, see |28j . 

Significant simplification follows for the case of a symmetric unit cell where the 
Fourier expansion ( l46l) i reduces to a cosine expansion with real valued coefficients. In 

consequence the matrix N is real and hence the basic conclusions obtained above for the 
case of a uniform half-space can be extended to the present case in hand. In particular 
by analogy with eq. (123]) it follows that the dispersion equation on the subsonic surface 
wave speed v = v s can be taken in any of the following equivalent real-valued forms 

det Z = <=> det (I + iS) = det B = tr adjB = 0), 

where adjB is the adjugate matrix. Note that the double zero of det B at v = v s makes 
the two other forms of the dispersion equation more appealing for numerical evaluation 

Of Vg. 

5.2.2 Examples 

Numerical results are presented for a laterally periodic half-space composed of layers of 
two alternating materials. Interfaces between the layers are normal to the surface of 
the half-space. Note that a periodically bilayered structure infinite along the periodicity 
axis can be seen as a periodically tri-layered with a symmetric unit cell and hence with 
a real matrix N, see §5.2.1. We assume bimaterial structures of equidistant layers of any 
two of the three isotropic solids: copper (Cu, Young's modulus E = 115 GPa, Poisson's 
ratio v = 0.355, p = 8920 kg/m 3 ), aluminum (Al, 69 GPa, 0.334, 2700 kg/m 3 ) or steel 
(St, 203 GPa, 0.29, 7850 kg/m 3 ). The surface wave propagation direction m, that is 
the direction of the wavenumber vector k = km (see (I4"T|) ). is measured by the angle 
ip between m and the plane of interlayer interfaces, so that if) = corresponds to m 
in the layering direction, i.e., normal to the interfaces. The values of wavenumber k 
considered are restricted to the first Brillouin zone defined by the cell length in the 
layering direction, which is taken as unity. 

Figure [1] demonstrates the computed surface wave speed function of wave 

number k for seven distinct directions of m. The curves displayed are azimuthal cross- 
sections of the subsonic dispersion surface v s = u(km)/k = v s (k,ip). The numerical 
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results show that surface wave speed decreases as a function of k and increases with i/j 
for all bimaterials considered. The dependence of v s on k at fixed ip is greatest for waves 
traveling across the interfaces (ip = 0) and least for waves traveling along the interfaces 
(z/> = 90°). 

6 Conclusions 

Use of the matrix sign function provides a new and broader perspective of the surface 
wave problem in anisotropic media. It straightens out the methodology of the underlying 
matrix formalism and offers a direct method to compute the matrices involved. Starting 
from the defining property of the matrix sign function ffTTj) we have obtained known rela- 
tions for the Barnett-Lothe matrices, the impedance matrix and the dispersion equation 
for the surface wave speed, see eqs. (TPS]) . fl2Tl) and fl23|) . These expressions are achieved 
using only the sign function of i times the Stroh matrix N = KN T K without specifying 
its method of evaluation. An integral representation f|T3|) for the matrix sign function, 
combined with the explicit structure of the Stroh matrix f )27|) . leads immediately to 
the Barnett-Lothe integral relations f )25|) . In this paper we have concentrated on using 
the matrix sign function for direct formulation of the dispersion equation, without dis- 
cussing further properties of the Barnett-Lothe matrices which underlie the existence 
and uniqueness considerations in the surface wave theory [2, H, 0]. We have constructed 
an explicit solution of the matrix Riccati equation for the impedance Z by rewriting 
the Riccati equation in a form (134"1) that involves the matrix sign function. Apart from 
providing for the first time a direct solution of the Riccati equation, the use of the ma- 
trix sign function shows how this nonlinear equation is intimately related with the Stroh 
matrix. 

Perhaps the greatest advantage gained by using the matrix sign function is that 
it provides a natural formalism for framing the problem of subsonic surface waves in 
laterally periodic half-spaces which are inhomogeneous along the surface and uniform 
in the depth direction. We have shown how much of the structure for the homogeneous 
case carries over to the case of laterally periodic materials. For instance, the conditions 
for surface waves (J63l) and the form of the block matrices in the matrix sign function 
( )62|) mirror their counterparts for the homogeneous case, eqs. (122|) and (125|) . respectively. 
Naturally, there are major differences between the problems. Conditions for the exis- 
tence of surface waves in periodically inhomogeneous materials have only been recently 
established [30] but the methods that have been proposed for finding them are not as 
straightforward as for the homogeneous case. The approach that we have presented for 
the laterally periodic case, being linked via the matrix sign function to the well known 
formalism for the homogeneous half-space, offers, we believe, a clear and logical route 
for finding surface waves. Future work will examine this calculation method and the 
properties of the subsonic waves in more detail and will also consider supersonic solu- 
tions in the upper stopbands and passbands - the extension motivated by the classical 
paper by Alshits and Lothe (3l| . 




Figure 1: The surface wave speed v s for three laterally periodic bimaterials made of equal 
volume fractions of any two of copper, steel or aluminum (Cu, St, Al). The propagation 
direction m varies from ip = (normal to the layers) to ip = 90° in increments of 15°. 
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Appendix: Sign function and related matrix functions 

The matrix sign function of M is closely related to other standard matrix functions. 
The matrix projector functions P + and P_ are defined 

P ±£=jo where C± = : M£ = A£, Re A <0, }, (66) 



and therefore 



signM = P + - P_, I = P + + P P± = -I± ^signM. (67) 



The projector functions can be expressed via integration, e.g. [32 



P ± = -L / dz(zlTM)- 1 (68) 
2tt« Jc± 

where C± counter-clockwise encloses the finite right (left)-half plane {z : Re (±2) > 
0, \z\ < R}, R arbitrarily large. The Schwartz-Christoffel transformation z — > (z — 
l)(z + 1) _1 reduces this to an integral around the unit circle 



P± 277/ 



— ! dz(zI-F ±x ) \ F = (M-I)(M + I)- 1 . (69) 

7™ J 



|*|=l-0 



While F ±l are singular if M possesses eigenvalues =pl, we note that (1681) is unchanged 
for M -»■ aM + i&I, a, b e R, a > 0, and flESJ) is invariant for F ->■ (M-aI)(M + a*I) _1 , 
a e C, Re a > 0, and hence (J69l) can always be made regular by selecting a in the right 
half-plane. The integral representation (fTBl for the matrix sign function may be easily 
obtained from eqs. (1571) and (|69|) . 

The disk function [lit] is defined such that if M = U^ 1 ^\ U _1 where the eigen- 
values of Ji, J 2 have magnitudes < 1, > 1, respectively, then diskM = U^* °ju _1 . It 
follows that 

sign M = I - 2 disk ((M - I)- X (M + I)) . (70) 
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